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Abstract 

1 We present a list of formulas describing classical radiation of the rank s tensor field from an accelerated 

point-like source in flat space-time of arbitrary even dimension d. This allows straightforward evaluating 
the total intensity and radiated momentum for any s and d algorithmically, by hands or with the help of a 
computer (e.g. with an attached MAPLE program). Practical application of formulas is limited, because, for 
s > 1, the energy- momentum tensor for the point-like source is not conserved. This usually means that one 
cannot neglect contributions to radiation from tensions of the forces that cause acceleration of the source. 

> 

Q\ | Traditionally radiation processes were considered as subjects of direct physical application and therefore 

were deeply investigated only in no more than d — 4 space-time dimensions, see [1]. Even for d = 2 and 3, 
where obvious applications exist, say, to sound waves in media, the theory remains badly represented in the 
literature. Only recently, after the strings-inspired multi-dimensional models [2] attracted increasing attention 
[3], some papers on multi-dimensional radiation began to emerge [4, 5, 6, 8, 7]. Of course, they are still too 
few to cover the field at the level of exhaustiveness, typical for the literature on Ad radiation. In this paper we 
make a step, which we think is necessary for a systematic study of physical effects. Namely, we provide a list of 
generic formulas, describing classical radiation for arbitrary dimension d and rank s of the radiated fields, what 
should help to reveal the underlying physical and mathematical structures. In particular, one can immediately 
read off the radiation damping force in higher dimensions from our results, for example by the method of [4]. 

To illustrate the main results of this paper, we use several Tables. We begin with Table I, summarizing 
the logic of radiation calculus. The left column is our main line, the right column enumerates related subjects 
which deserve deeper physical consideration and are left for a detailed presentation of radiation physics in a 
more sophisticated review article. 

According to Table I, in order to find radiation intensity one needs to perform the following chain of 
calculations: 

• Solve wave equation for a point-like source of the rank-s field, moving along a world line z m (t), 

□Au-m.O*) = ju^.-.u^^x- z{r))dr, (1) 

and pick up the slowest-falling contribution at large distances. For even d it is given by a simple formula for 
the retarded Lienard-Wiechert potential: 



X 



where = d T z^ is source's d-velocity, u 2 = 1, evaluated at the radiation-time moment t 1 = z (t) which is 
determined by the condition R 2 = 0, where i? M is the d- vector with components W = x^ — z^(t). Here r is 
source's proper time. In particular, 

— = (3) 
dx» (Ru) ' 
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We also introduce the notation n M = (1, —ft) = 5£-, where R is length of the spatial part of R^, R = (R 1 ) 2 
(n M is not a d- vector in contrast to we define it in the laboratory frame). 
It is convenient to re-write (2) in condensed notation: 

A ^ = ( jUY 1 jj = ja-5 E (d q T S)A^ and d,A rad = E (d^)A^ (4) 

^ ' q=0 q=0 

where U = (nu), S = u fil . . . u^ s , p = ^p, and find A v q from recurrent relations 

< 1 = i(a r ^ + ^_ 1 ) (5) 

The first few A\ are listed in Table II. 



Check transversality: in the leading order in 1/R, when derivatives of R are neglected, 



d-2 



^ M ^-m. = ((^)^J 2 (6) 

does not vanish for s > 2, and this is not cured by subtraction of traces. The physical reason for non- 
transversality is the disregard of radiation from tensions of the forces that cause acceleration of the source: only 
taking into account all the tensions makes radiation problem well-defined for s > 2. For s > 2 the formulas of 
the present paper provide only a part of the full answer. 

• Build up the energy-momentum tensor in the wave zone. It is equal to 

« - { (jk ,% ) * " ^ {§ A ' (7 > 

In applications, certain linear combinations of such stress tensors arise. For instance, for the scalar waves 

rp{scalar) 7^(0) (o\ 

This is because only the spatial components of all non-zero spin fields have any physical meaning (e.g. survive 
in physical gauges), thus giving rise to the overall minus sign of the kinetic part of the energy- momentum tensor 
as compared with the scalar case. For gravitational waves 

T (grav) _ rp(2) _ 1^(0) (q\ 

Similar redefinitions are also needed for higher spins s > 2. These linear transformations can be easily made 
using the results for T$ read off from Table V. 

• Evaluate the radiated momentum flux through the sphere of radius R. It is equal to 

dV^ s = -j> {T^nrfx } dS (10) 

where dS is an infinitesimal element of this sphere (i.e. integration runs over the sphere). Then, one finally 
obtains (using dx° = Udr = (nu)dr) 



dr 



|s = - <L T^n t R d - 2 Udn d _ 2 = 



= £ («*...«"•)) f^^^^H^^pf^,) (ii) 

q',q"=0 k=0 

where dttd-2 is the solid angle in the d — 1 space and are certain functions of various scalar products 
(d l T u v d™u u ) with I + m < d — 4 — k. In particular, the radiated energy loss is (dt' = u°dr = jdr) 

d-4 

= - E (d k T j)Pt dls) (K) f ° r C ° nStant ^ (1 -" 2)_1/2 P^\ K ), (12) 

^ fe=0 
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We warn once again that due to effects like (6) and (9)- (8), this quantity represents only a part of the full 
answer for s > 2. Moreover, as follows from Table V, it can be even negative (!): for example, for fields with 

asymptotically large rank s the leading contribution to the radiated momentum, d ^ 3 behaves as — (su 2 ) 3 _1 ~ 

(— )i , provided u 2 < 0. For s = 0, this is cured by changing the common sign, see (8). However, for s > 3 the 
procedure has to be more sophisticated. 

Sometime also contractions I d \ s — are considered [6], but they do not have any direct physical 

dV^ dV° 

meaning unless d ^ s ~ u^, as happens for d = 4. In this latter case, I d \ s = d p s . 

• Calculate angular integrals over isotropic (light-like) unit vectors in (11) (although is not a d-vector, 
these integrals are Lorentz-covariant, since (11) is Lorentz-covariant)), 



/ 



dn d . 2 (n) ~Pr«,(u" 1 ...u"™) (13) 



They are expressed through the trace-eliminating projector Pr<j (since n 2 = and contraction of any pair of 
^-indices at the l.h.s. should give zero): 



spin s = 2 
spin s = 3 
spin s = 4 



l 

d+A 



Pr d (u AJ M, y ) = u^u v - ^ri^ 

Pr c ;(ii^ 1 ii^ 2 ) — u^ 1 u^ 12 u^ L3 u^ l4 

Vtiifj.2 u fi3 u Ln ~b 5 permutations^ + ^pjy^p^ {vtn^V^tn "b V^i^V^^a "b ^i ^4^2/^3^ 



For generic m we have: 

/ ^2fc permutations 

[ ^C 21 fc (rf + 2m-4-2fc)!! ' 7 " 

"r<j(u Ml . . . u Mm J — ' (d + 2m — 4)!! m mi ■ • ■ "^-it ^I^^Hi^^t+i • • • ^ 



Y {2k — 1)!! permutations 



This projector can be conveniently re-written through generating function: after contraction with m copies of 
a d- vector x we get 

Pr d ((x«r) = ER^fc - l)UC? k {d + ^ 4 _~ 2 * )!! { X uT- 2 \x 2 f = 

k=0 

_ y fc (d + 2m - 4 - 2fc)H m!(^) m - 2fc (x 2 ) fc 

(d + 2m-4)!! 2 fe /c!(m-2fc)! ( ' 

The sum of coefficients in this formula is needed to define the normalization factor in (13): 

« . (d + 2m-4-2k)« (d+2[m/2]-3)n(d+2(m-[m/2])-4)!! 

c(m) = g h'cbp* - i)aL+__L _ i _± j L__ L m 

Here [m/2] means integer part of a number, e.g. [3/2] = 1 and 3 — [3/2] = 2. 

Rewritten in terms of the generating function and with normalization factor restored, the angular integral, 
(13) is (see Table III): 

(nx) m dn d _ 2 (n) _ S n _ 2 / x _ 
{nu) d+ m -2 Cd (s) d V ' J 

2d/2 ^ y\-y (d + lm-A-m m! , )m _ 2j , 2)j (ljJ) 



(d + 2[m/2]-3)!! ^ v ; {d + 2m - 2 [m/2] - 4)!! ?Jj\(m - 2j)\ 



d — 2 

^k„ 



It remains to substitute the vector x by x = txo + X)fe=o a k9^u and pick up the coefficient in front of the 
relevant combination of t and a*,. Note that this calculation uses the generation function which provides an 
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additional combinatorial coefficient that has to be taken into account manifestly. E.g., in order to get the correct 
coefficients in front of the term d^u . . . d^ n u with different values of k\, . . . , k n , one has to divide the generating 
function by nl. 

The final results can be equally expressed via scalar products Uij = d l T u v d" l u v and through the Frenet 
curvatures K m and their derivatives d l T K m . Hence, the last step: 

• Re-express Uj through Frenet curvatures [10] which parameterize the moving orthonormal basis, associated 
with the world line x M (r). It is formed by d- vectors /j, = 0, . . . , d — 1, 

N^N") = rT (17) 

with — u (i.e. = u^) and other vectors constructed recursively from the condition 

From r-derivative of orthonormality condition, 

(d T N^NM + NM (drN^ = . 

It follows that at v < fi 

/V = (d T W)NM = -N^(d T N^ = K^-nf 1 

and 

d T N^ = -k^ +1 N^ + V + k^N^ 1 (18) 

Parameters fc M (they are not d- vectors!) depend on the shape of the world line x^(t) in the infinitesimal vicinity 
of its point and are called Frenet curvatures. Expressions of scalar products Uij through these curvatures are 
given in Table IV. 



Results 



All these steps are easily made with the help of MAPLE or Mathematica, see Appendix at the end of this paper 
for a MAPLE program, [9] . As an application of the algorithm we obtain the explicit formulas for the radiated 
d-momentum dV^ s as well as somewhat simpler formulas for P^ s \ the radiated intensity at constant 7, (i.e. 

2 2 dV^ 

K 2 = tzi =const, Ki — 0, i > 2) see (12), and for the contractions Id\ s = J^" , [6] for "realistic" values of 
d = 4, 6, 8, 10 and generic rank s. First, we list results for the case of d = 4, 6, 8, when the formulas are relatively 
simple, and one can understand their general structure. Much more involved formulas in the most interesting 
case of d = 10 are put into a separate section. 
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4- 12s 



Cut 



dr 



= TT 



19 



- (2s - 3) 5 



8tt 



16tt 



16^ 2 



w 4 ^ + ^-(1 - 5s)il 2 u^ + -^-(2 - ls)(uu)u^ + ^r(7s - 4)u 2 u» 



15 



35 



105 



dr 



16tt 3 /7s 3 135 2 838 \ , 32tt 3 
— (lO "14^ + 35- 5 ) U + — 



2 101 41 . . , 



16tt 3 



s 2 59 20 



+ — { 'J + i S - i J " V + ^ 15 - Y + 189 ^ 15 



16tt 3 



f 32tt 3 / 6s 2 19 283\ ....... 64tt 3 / 1 s\ .... 

{ — + T S - 297 j + — (27 - 5 j tttt > ^ + 



+ 



64tt 3 /2s 2 



21 



65 , 
2s+ 99 |/f 



64tt 3 /l 2s 



21 V 9 5 



uu > + 



32tt 



17 



. s — — 1 (mi)it 
35 V 27 M ; 



(d|s) 



( 4| S ) _ 12a -4 2 

~0 ~ o 7TK l 



p(6|s) _ 

mi — 



V I 4s^ - yS 



28 32 



15 



8tt 



4^ + T ^(5s-l)« 2 « 2 



^ 8|S) = lk {7s + ^ + fi^ 49s3 - 549s2 + 10045 - 216 ) K i - Wn3 ( y - a + ) ' 



8tt 3 



64 



12s - 4 



'4|s 



/„: = -TT 2 ( 4s 2 -™ s + 2 ^\k{ + y^(5s - 1)(k 2 + 



'8|s 



16tt 3 
105 



49 , 167 n 406 276 



172 



I S T + T S_ ll K i _ I 21s - 97s + -g- J K\Ki - I 77s — 180s H — — ) 



109 



2(7s 2 - 9s + 4)kJ«i + (7s - 1) {^i^k 2 + _ «i) 2 + ( 2 «2Ki + kik 2 ) 2 } 
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Results for the most interesting case of d = 10 



dV" ms = f 4tt 4 
dr \ 315 



6466 , 236237 2 1609666 9800 \ . 8 32tt 4 



19s 4 - — — ^ + 



33 



99 



7 



200 2 2147 6230 N 

Vs" s A s 
9 33 429 J 



+ 



Wt, ( 3 61^ 2 8316^ t 11690 \ 1; . 2 



35 



27 297 1287 / 



64tt 4 / , 1037 9 32653 4865\ ....... 

" ^ : < " ^ + "IT I s " ™ + " 429 J " M ~ 



1671 " 4 An 2 611 512\ „ 4 32tt 4 /17 2 irl 34 \ . . .... 2 64tt 4 / 2 53 100 \ , 

-105" ( 19S - ~ S + 33 J U ~ "IS- U S - 1515 + 21 ) {UU) - K)5- I 5 " Y S + 33 ) MW 



256tt 4 / 2 263 119 \ 2/ 32tt 4 / 2 134 25 \ 2 -2 64tt 4 / 2 74 15\ 

~35" S T08 S + 198 j " W - T05" S " — S + IT U u + K)5- S " Y S + IT " + 



64tt 4 ^ 2 ,„ ,„ w ... w .::, 32tt 4 /l 
+ ^ (9S 495 + 10)(UU)(UM) + T05" 1 9 



s u 



}"-{ 



32vr 4 /8s 3 593 2 , 72173 ^ 5670 \ 



63 V 3 n " S + 429 " 143 J 



u 4 (uu) + 



647T 4 , 



256tt 4 / 2 133 4 5 3 \ 128tt 4 / 2 109 147\ ....... 



, 829 2664\ . . 2f . .... , r _, 

11 s 143 / 



320tt 4 /_ 208 \ , i2( ,,:. A , 64tt 4 ^ 3 ^ f ;.^\ ^ ( / 128tt 4 f 3 785 ^ , 30843 ^ 2065 \ . 6 



819 



33 



+ 1 - — s u 2 (iiii) + — — 2s - — (ilu) } + 



189 
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\ 189 44 S ' 572 " 143 J 



320tt 4 / 2 53s 216 \ , N? 640tt 4 / 2 1569 \ 2/ - N 64tt 4 / 2 I 23 ", 

+ ^8?r S " IT + 143 M - "W (" - 45 + 1430 ) U ^ ~ W 4 * + * + 143 1 U " + 



+ - 



128tt 4 / 2 



189 V 11 

32tt 4 



945 





f 256tt 4 




\ 4158 


141 


1020\ . 






TT S4 


- 143 r 



™ 2 795\ ....... 64tt 4 / 3 

22s 2 - 167s + - u 2 (mi) + — ( — - s ) («u) } u 



945 V 44 
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32tt 4 
315 



19 4 2747 , 152021 2 369827 1 \ 8 /9 , 2 „ 00 „ \ 
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167T 4 

945 



s ,^no 2 53244 133048\ B 2 ,„„ r 2 106818 300280\ 
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4,%2 
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- ^99s 2 - 931s + k\k\ + ^36s 2 - 220s + k\k 2 k 2 + ^324s 2 - 1780s + k^k 2 - 

- (^468s 2 - 1160s + Kl «!« 2 + (l44s 2 - 568s + ^) k 2 /^ - (l71s 2 - 413s + ^) A*- 
2 .oo 866 \ 422 A 7 4 2507 3 125989 2 3644975 111744\ 8 



54s^ - 488s + ^ j ^ - ^-j-* - + -77-^ " + 1ST 

- (l08s 3 - 1506s 2 + ^g?) - (W - 822s + ^ ( 43 2s 2 - 2752s + «? 

- (W - 268s + ^) K ^ 2 - (450s 2 - 3714s + ^) k 2 k 2 k 2 + 
+2(9s — l)[(«i — 3ki«2 — 3kik 2 K2) 2 + (3k!K 2 + ik\k 2 + Ki'k 2 ~ kik 2 ) 2 + (2nik 2 n 3 + KiK 2 k 3 + 3K1K2K3) 2 

— k\k\k\{2k% + Kg + K 2 ) — 2kiK 2 K 2 (3kiK 2 + Kl«2 + 3Ki«2)]] 



KlK 2 K 2 — 



Note that, in the expressions for the radiated momentum the coefficients in front of terms of degree k in U%j 
are polynomials of degree k in rank s. Relative simplicity of the coefficients of these polynomials implies that 
a general formula may exist for them at any dimension d and for arbitrary rank s. However, we stress once 
again that these expressions for higher ranks are of restricted use, because of necessity to take into account 
contributions from other energy-momentum tensors. The second thing one has to do in order to separate the 
contribution of the given spin s from the expression for the rank s is to add proper combinations of energy- 
momentum tensors for lower rank fields. This is done straightforward using formulas we obtained here. 

The examples of these formulas, with illustrative purposes, are collected in Table V for low values of 
s = 0, 1, 2, 3 for dimensions d = 4, 6, 8 when the formulas are more comprehensible. Much more involved (but 
interesting) case of d = 10 is listed after that table. The Appendix contains an example of the MAPLE program 
for calculating various quantities obtained in this paper. 
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Tables 



TABLE I 



Radiation: calculus and problems 



accelerating point-like body (source) 
I 

radiation field 

I 

even dimensions 



\ 



separation of radiation from co-moving field, wave zone 



strong violation of Huygens principle at finite distances 
in odd dimensions 



d-4 

retarded field = (^0 T ) T «i — difference between observation and radiation times 



field tension 



transversality violation for s > 1 

d-2 

other sources of radiation can not be neglected 



stress tensor 



« ) {(RU) 9 !-) 



d-2 s 2 

2 u"i...u"« 



(Ru) 



i 

full intensity 
I 

polarization, angular and frequency distributions 



picking the definite-spin contribution, 
radiation from other sources to cure T M1/ non-consevation 

radiation friction 
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TABLE II : A\ 


9 = 


9 = 1 


q = 2 


q = 3 


9 = 4 


field in d = 4 

p = 


i 

u 


- 


- 


- 


- 


field in d = 6 

p=l 
tension in d = 4 


u 

U 3 


1 

U 2 


- 


- 


- 


field in d = 8 
p=2 
tension in d = 6 


-UU+3U 2 

Tt*\ 


_3U 

u 4 


U 


- 


- 


field in d = 10 

p = 3 
tension in d — 8 


-UU 2 +WUUU-15U 3 
V 


-4.UU+15U 2 


6U 

U 5 


1 

U 4 




field in d = 12 

p = 4 
tension in d = 10 


~UU 3 + 15UUU 2 + 10U 2 U 2 — 105UU 2 U+105U 4 


-5UU 2 +60UUU-W5U 3 


-WUU+45U 2 


wu 
1W 


l 

m 
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TABLE III 



Formulas (16) 



for angular integrals 



d = A 
(p = 0) 



d = 6 
(p = l) 



d = 8 
(p = 2) 



d= 10 
(p = 3) 



m = 1 



A-k(ux) 



^-(ux) 



iffl(«x) 



^-(«x) 



m = 2 



f|Li(lO(«x) 2 -^ 



m = 3 



^^6(nx) 3 -3(W):r 2 ) 



8(«:r) 3 - 3(ua;)x 2 



lO(urr) 3 - 3(ux)s' 



12(ua;) 3 - 3(W):r' 



m = 4 



05 ' 8 ( ux ) 

-8-6(ux) 2 x 2 + 3(x 2 ) 2 



-10 ■ 6(ux) 2 x 2 + 3(x 2 ) 2 ^j 



||g±(l4-12(«:r) 4 - 
12-e(ux) 2 x 2 + 3(x 2 ) 2 ^j 



m = 5 



105 



■ (l2 ■ 10(mx) 5 - 
-10 • 10(ua;) 3 x 2 + 5 ■ 3(ua;)(a; 2 ) 2 ^ 



945 



(l4 ■ 12(ux) 5 - 
-12 • 10(nx) 3 a; 2 + 5 ■ 3(ux)(x 2 ) 2 



iW ( 16 -14(«x) 8 - 
-14 ■ 10(nx) 3 x 2 + 5 ■ 3(nx)(x 2 ) 2 ^ 



m = 6 



■ (l6 ■ 14 ■ 12(mx) 6 



-14 ■ 12 ■ 15(ux) 4 x 2 + 
+12 ■ 15 ■ 3(mx) 2 (x 2 ) 2 - 15(x 2 ) 3 



jUfgg |^18 ■ 16 • 14(ux) 6 - 
-16 • 14 ■ 15(ux) 4 x 2 + 
+14 ■ 15 ■ 3(mx) 2 (x 2 ) 2 - 15(a; 2 ) 3 



m = 7 



^(l8.16.14(ux)*- 
-16 ■ 14 ■ 2l(uxfx 2 + 
+14 ■ 35 ■ 3(nx) 3 (a; 2 ) 2 - 15(«:r)(:r 2 ) 3 



^ (20-18- 16(ux)?- 
-18 • 16 ■ 21(ux) 5 x 2 + 
+16 ■ 35 ■ 3(nx) 3 (x 2 ) 2 - 15(ux)(x 2 )^ 



m = 8 



(22- 20-18 .16(«x) 8 - 

-20 ■ 18 ■ 16 ■ 28(nx) 6 a; 2 + 
+18 -16-70 ■ 3(«x) 4 (a; 2 ) 2 - 

-16 - 28 ■ 15(«:r) 2 (:r 2 ) 3 + 105(:r 2 ) 4 ) 



m = 9 



^f%^ (24- 22- 20 .ISCua:) 9 

-22 ■ 20 ■ 18 - 36(nx) 7 x 2 + 
+20 ■ 18 • 126 - 3(uz) 5 (a; 2 ) 2 - 
-18-84-15(«a;) 3 (x 2 ) 3 + 



+9 



105(mx)(x 2 ) 4 ^) 
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TABLE IV : 


Scalar products through Frenet curvatures 


Unit vectors 7V», nMnW = rf v 


order 


u 2 = 1 


JV(°) = u, 


u(o) 

A" = KiA^f 1 ) 


order 1 


uii = 
u 2 = -n\ 


• 


U) = Kl Af(°) - K2 Af( 2 ) 


order 2 


uil = 4 
uil = —K\ki 

"2 4 2 2 • 2 
U — ^\ ^1^2 ^ 1 


KlK 2 \ Kl 1 / 

1.(2) 

TV = K 2 iV (1) - ^ (3) 


order 3 


uii = 3/ti/ti 
itii = —K\k\ — k\ + 4 K 2 
ilil — —ki'ki + 2k^ki — kiKjKi — 4 K ik2 

if = — k\ - 2n\ki + 2ki4^i + 9k 2 k 2 — ^n\k\- 

4.K l ^2 ^1 ^2 ^1^2 ^1 2^^/^2 ^1^2 ^1^2^3 


iv< 3 )- 1 (k { 2k ^ + k -)ti+ 

K1K2K3 ^ \ K\ K 2 > 

+ (4 - k\ + ^4 + ^ - '*A'u + 4(^ - ^-)u) 

1.(3) 

AT = K3 AT( 2 ) - K4 Af( 4 ) 


order 4 


uii = 4kiKi + 3k 2 + k| — k 2 ^2 
itii = — 6Kf/ii + 3ki/v|ki + 3k 2 K2«2 — ^l^i 
iiii = — k\'k\ + 4/tf/ii — 3ki4ki — k 2 K2«2 — 

3Av^Av-^ — t - 3^2^1 — t - ^1 2^^^2 — t - ^r^2 ^1^2^3 
-iiii = — itiKi + 6K 2 /ilKl — KiK?,K 2 ki — K 2 K2K 2 k2 — 

— 4 K2K3K3 — 3K2K1K1 + 9k±4 + 444^1 + 2nfK2K2 — 

— K\K\k\ — 2k 2 K2K2 + K\K^k\ — 3«fKl — K^k\ — 
—2K2k\K\k2 — 6K2«f «2 — 3KiKiK 2 — 
7/ — ~~ ~ 9 t\j 2 Kg ^1 ~ — ^1 ^2^3 ^4 — 6k 2 ^ 1 Ki ^2 ^2 
— 1444^1 — — 4k 2 K3K2K2«3— 

— 6kiK2KiK2K§ + 30k 2 k|k 2 — 4 K 2 K i — 2n\n\n\ — 
— 9k|k 2 + 30Kf«i«;2K2 — \2K^k\K\k2 + 8«;f «i+ 

+ 16k 2 K 2 — 9k 2 K 2 — K 2 K2 + 6K1K2K2KI — I8K1K2K2K1 — 
— 6kiK2Klk2 + 6KiK2 K 3^1 + 2k 2 K2K 2 K2 — 6K2K1K1K2 + 

+24kiKi« 2 — 12k 2 kiKi — 2«f K2K2 + 6«ik|ki + 

. 2 

+ 2k 2 K2K2 — k 1 


A^W = 1 [ iT+ fs^ 1 +2« 2 ^0 ti+ 

\ ( J2 ^2 _2 1 0K1 fi K 2 k 2 
+ ^1 - «2 - «3 + 3— - 0^ - 2^- 

_^M2 _ 2 K 1 K 3 _ K2K3 ] 1 ( «1 I g K l I 

KiK 2 K1K3 K2K3 y yKi ' Kf 

+( K 2 + 4 + 4)^- (24 + 4)%- («? - 

. 2 • • 2 

g KjKi _ KjK 2 _ 2 ^1^2 _j_ «l«2 _|_ 2^1 ^2 

k\ K±K2 K\K2 K1K3 K\K2 Kl«2 
2 K|K3 + K1K2K3A ^ +k 2 /Kl K 2+fi .2 + 
K^K3 K1K2K3 / J- V Kl K 2 O 

2 K 2 _ O K1K2 _ K1K3 1 K2K3 j ^ 

kik 2 K1K3 K2K3 y 

1.(4) 

A" = K 4 iV (3) - /t 5 AT (5) 
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TABLE V: 





Radiated momentum flux 


Particular Projections 


1 el_LlJ\. U LL.1I0UI 

vector 
rank 2 tensor 

gravity 
rank 3 tensor 


' *|o - 4tt - 2 p 
dr — 3 U U 

dT>^ 

dr — 3 " u 
d"P^ 

aP 4|2 _ 207T - 2 u 

dr _ 3 " u 

d"P M d"P^ 

aP 4| 3 _ ?M l a '4|0_ 227r - 2 7 ,u 
dr ~~ ; 4|2 2 dr ~~ 3 " U 


P (4|0) - Lin - -^-K 2 

p( 4 |l) _ r , _ 87T..2 

p(4|2) _ j _ 20tt 2 

p(4|ff) _ r _ r _ 2271-2 

p(4|3) _ x _ 3271-2 
M) - 7 4|3 - —«1 


rank s tensor 




P (4|s) - 7\,, - 12s -4 7rK 2 
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Radiated momentum flux 


Particular Projections 


dV / 10 - 8 f + *g u 2 u» + 32 f {uu)u» ^ u 2 il^ 


7 6 | = f^(20Kf + 7K 2 + 7^ 2 ) 


^ 1 = 32x 2 -2^ ^tyw + 16V tf-p 


d(6|0) 87r 2 /,i 4 , 2 2\ 

7 6|1 = ^(9k? + 7k 2 + 7k 2 k 2 ) 


-jili = l^-u^ - ^-u 2 u^ - ^-{uu)u^ + ^-u 2 u^ 


P («W = Ifgi^f + 2k 2 k 2 ) 
J 6|2 = ( - 13k$ + 63k 2 + 63k 2 k 2 J 


dV I 19 = 2 °f « 4 «" "if ™f (««)«" + 35^ 


P (6|2) = ^ + 9. 2 k 2 ) 

4|g = 4|2 - \h\0 = 

= w ( ~ 46k? + 119k? + 119k 2 k 2 ) 

105 y i i 1 / 




P (%) = ^( 6 «* + l9^) 
J 6|3 = TM ( - !67k| + 98k? + 98k 2 k 2 ) 

p (6|3) = _^ (19k 4_ 14k 2 k 2 ) 


6|s ^2 r 19 /o„ Q^2^ „\4„,u i 8w (-i r„\„",2„,u , 
rfT = 7T l-g (Zs — 6) J U H j-g- 1 — osju + 


r ^2 / a „2 124 „ | 32 \ ,A , 


+±§^(2 - 7s)(«ti)t*'* + ^(7s - 4)ii 2 u^ 


+ ^(5.s-1)« 2 + ^(5 S -1)k 2 « 2 ) 

P (6|s) = -7T 2 (4s 2 - f s + ff) 4k| + 
+ ^(5s-1)k 2 k 2 ) 
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d = 8 



dT"* 

a/ 8|0 



dr 



8|1 



_ / 1607T 3 -6 I 12448^ 3 / • ■ ^SSSir 3 .9- .9 AKfiAir 3 .9/ . -\ 39ir 3 -21 

— \ 9 U H 945 



544tt 3 

|_ 2079 " v / 189 j - 1 |_ 2079 " 1 189 J " 945 - 

= {-I^ + »i(^) 2 + ^u 2 « 2 - ^u 2 («ii) - ^ 

, f 782727T 3 - 2/ ■ -;\ _ 14087T 3 • ai _ f 29888tt 3 -4 , 832tt 3 

r \ 10395 945 UU J " I 10395 ' 



al a 



ao27r° -2"2 
189 " U 



_ f 6952tt 3 - 6 I 14576tt 3 / ■ ■■\2 , 2432; 
dr — \ 315 " ' 945 y aa > ' 189 

dy, 8|3 _ / 5927^-6 _ 100167T 3 / ■ -;N2 
dr ~~ \ 315 U 945 [ Ua > 



>(8|0) 



-t n 



= -^(3K 2 K 2 K 2 + 108 K f+ 
+64k|k| + 3k 2 k|) 



+314^k 2 + 39k 2 k|) 



j(8|fl) _ 72tt 3 „2„2„2 , 256tt 3 
— ~35~ re l K 2 re 3 "T 105 h 

I 5536tt 3 4„2 , 72tt 3 2„4 
315 - K 1 K 2 "r _ 35 _K 1"'2 



= ig^(60^ 2 K 2 -411^ 
+314/tf«; 2 . + 60k 2 k|) 



dPi 



dr 



J" 16tt 3 / 7 

I 9 l ] 

167T 3 /_ 
3 I 



- ¥ * 2 + If- * - 5) - 6 + ^ H 2 + # - - i) (-) 2 + 

+ 1^-11 « 2 « 2 + 32rr 3 (fi - f + ^) « 2 (^)+ 



,(8|«) 



167T J 

105 



(7s - 1)(k 2 k|k| + K 2 /t|)+ 



+ lffi(i- S )^ 2 }^ + {^(-^ + f S -if)^) + ^(i-|)^} 
+ {^( 2 #-2 S+ i)^ + ^(I-f)^}^ + ^( S -lI) (uu)il- 



+|fg (49s 3 - 549s 2 + 1004s - 216)k?- 



16tt 3 
1155 



276k? + ^^k^k^ + ^^-n\k\ - 88k?ki + 11 + - k^ 2 + (2k 2 «i + «ik 2 ) 2 } 



32tt 3 
385 



_ . ft 2816^9 ^-^-^^2«2 iif222 / 2 "\2 • •\2 s l 

"4«i H 7p^~ K l K 2 ^ 7~ K l K l 7T K l K l + ll \ K l K 2 K 3 + ( K 1 K 2 — K l) + (2«2Ki + K1K2J ) 



h\2 = 



16tt 3 
H55 



-254k^ + k^k 2 , + ^-n\k\ - 132kJki + 143 {k 2 K2 k 3 + ( K i K 2 - ^i) 2 + (2k2«i + K1K2) 2 } 



z 8|g 



8rr 3 
H55 



-232«f + l9 ^ 88 K 4 K 2 + ^5 K 2^.2 _ 352k 3^ i + 297 {k 2 «2« 2 + (mk 2 . - ^i) 2 + (2k 2 «i + kik 2 ) 2 } 



h\3 = 



16tt 3 
1155 



-1651/tf + — fc^/t 2 - ^^k 2 k 2 + 220k 3 ki + 220 + - k^f + (2k 2 ki + K1K2) 2 } 



16tt 
105 



3 r 



49 



167 



406 



21s 2 - 97s 



172 



4 2 

K 1 K 2 



77s 2 - 180s + 



109 



2 ■ 2, 
^l^l~r 



+ 2(7s 2 - 9s + A)k\R\ + (7s - 1) {k 2 k 2 k 2 + (kik| - k x ) 2 + (2k 2 «i + «ik 2 ) 2 }] 
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d = 10 



<#Toio f 1120tt 4 1088tt 4 . . 8192tt 4 .. 4 32tt 4 ....2 28480tt 4 . ,, . ... 2 5344tt 4 , .., 

— ^ = <^ ?i 8 uti 2 m 4 H u H m 2 (mm) 2 H m 4 m 2 - 

dr \ 27 315 V ; 3465 945 429 1287 



8896tt 4 ,,. .... 2176tt 4 .. 2/ . .... 1280tt 4 , . ...... .... 160tt 4 . 2 -2 164tt 4 ....... 128tt 4 . . .... . .... . 

42CT " ^ UU ' 495"" ^ UU > 693~^ U ^ U ^ ~ ~m~ U U + ~1t~ U ^ + ~q^~ ( - uu ^ uu > < u + 

f647r 4 ,;.», 38656tt 4 . ..... .. 896tt 4 :: , 2880tt 4 ... ... 320tt 4 . .... 189447r 4 .. 2/ 

+ < m (uu)(uu) H u (uu) -\ u (mm) M (mm) m (mm) > M^— 

\ 693 v ; 9009 v A ; 1287 y ' 143 v ' 819 v ; 3003 v ; ' 

f377607r 4 „ 2560tt 4 x2 256tt 4 , ■■ 33472tt 4 2/ • 2624tt 4 2 ,1 „ 
- • ' ( uu ) + ~~^u7-(T '■ " " ' - (««) - ^7T" " «"+ 



t 3861 1001 v y 2079 v y 9009 v y 9009 

/ 64tt 4 . :: 33920tt 4 . 2 . .. 1 ... M f 1472tt 4 „ 2 1856tt 4 . ... 2176tt 



.4 



4 I "-M 



< MM H M (Mii) > il — < M MM M > U 

\ 693 9009 v '} \ 10395 10395 9009 



d ^w\i [ 2240tt 4 704tt 4 . 7552tt 4 , 256tt 4 ::2 15680tt 4 .,, . ... 2 320tt 4 , 2 

5440tt 4 63424tt 4 ,. .,. 640tt 4 , . ...... .... 7360tt 4 . 2 .,2 7424tt 4 . 2 , .. •• . 128tt 4 , . .... . -I „ 

+ 1? " W + 1Q395 " + — HW + 2Q79 " u - 2Q79 " (™) 2~j— (mm)(mm) I m^ 1 - 

f 1216tt 4 ......... 4352tt 4 . ..... .,. 1984tt 4 . 2 . ...... 50560tt 4 ... . ... 8000tt 4 . 2 ... .... 112064tt 4 .. 2 . 

- 1^079- (uu) + ^I7r (uu)(uu) - "TooT" {uu) - ^28T U {uu) + ^mT u {uu) + ^oo9- u {uu) ) u ' 

f592007r 4 B 3200tt 4 x2 128tt 4 , .. 58048tt 4 2/ . 53632tt 4 2 2 1 „ 

+ <^ M b (uu) (MM) H M (MM) M M ^ M M - 

\ 3861 819 1 ' 231 1 ; 9009 v ; 27027 J 

f 512tt 4 . .... 139 5 207T 4 . ...I ...^ f 128tt 4 „ 2 64tt 4 . ... 4288tt 4 . 4 1 

_ „^„^ MM — — — — — — M (MM) > M + < — — — M + ——MM + — — — M > M 

\ 2079 27027 v ' ) \ 2079 693 27027 J 



rf:P io|2 fl600007r 4 . 8 2048tt 4 , . ; .. 2 10016tt 4 , 544tt 4 -,.2 1429504tt 4 . 2 . . ... 2 1424480tt 4 ... 2 

^ = <^ M 8 H (MM) 2 H M 4 M M 2 (MM) 2 M 4 M 2 + 

dr \ 2079 315 V 7 1485 945 9009 1 ' 27027 

210304tt 4 . 4 ...... 40967T 4 .. 2 ,.. N 59776tt 4 , . ...... .... 74464tt 4 . 2 ., 2 70208tt 4 . 2 ... ..... 3328tt 4 .. .... ...... 1 „ 

+ ~m^ U {UU) + ^079"" {UU) + -M6T {UU){UU) + ~W39^ U U ~ ^039^" {uu) {uu) {uu) f U 

r 26247r 4 /; .:: N 371207T 4 89344tt 4 . 2 . . :: . 469120tt 4 . 4 . . ... 122560tt 4 . 2 ... .... 476800tt 4 ,. 2/ . ..J . „ 

_ { [UU) H ——— MM MM) — — M MM — — M (mm) H — — W MM) H — — — — — — M (mm) > M F + 

\ 2079 v ' 27027 v A ' 27027 v ' 9009 v ; 27027 v ; 27027 v ' ) 



f499847r 4 . R 1888007T 4 ,. ... 2 25607T . •• , 37952tt 4 . 2 . . ., , 57536tt 4 , ,1 ..„ 

+ { r~rzrz M b — — (mm) 2 — — (mm) + — — M 2 (mm) — — M M ^ M^- 

\ 2457 27027 v ; 2079 v ; 3861 v ; 9009 J 

f 1216tt 4 . •• 11392tt 4 . 2 , . ... "I : .n f 2752tt 4 „ 2 3776tt 4 . • 66368tt 4 .A 

— { MM M MM) > M + < U H MM H M > M 

\ 2079 1001 v ') \ 10395 10395 135135 J 
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dVw\ g f 2031207r 4 o 288tt 4 - n2 16480tt 4 4 16tt 4 -2 1728544tt 4 2 , x2 1480592tt 4 , 2 

- = <^ m s h utir H u u u 2 (uu) 2 u 4 u 2 + 

dr \ 2079 35 1 ; 2079 27 9009 v ; 27027 

250336tt 4 . 4 ...... 433287T 4 ., 20992tt 4 . . ...... ; ., 111527T 4 .,...2 74528tt 4 ....... 1216tt 4 . . .... . ..... . 

+ ^861-" {UU) + ^0395-" {UU) + ^155- (mi)(wi) + ^485"" u " ^0395^ " 105- (tm)(tttt) r "' 

f27207r 4 ......... 95104tt 4 98752tt 4 .,. . :: . 559840tt 4 . ,. . ... 127840tt 4 . 2 ... .... 562048tt 4 ,. 2/ . . , , 

_ S (uu) H — — — — (uu)(uu) — — - u m — — — u (uu) H „„„„„ it tra H — — — u (uu) > 

\ 2079 v ; 27027 v A 1 21021 K ' 9009 v ' 27027 v ; 27027 v 

757767r 4 . fi 223360tt 4 ... x2 128tt 4 ,.... s 315872tt 4 . 2/ . 8032tt 4 . 2 ,. 2 l ..„ 

u° (uu) (uu) H u (uu) u 11 > u^ 1 - 

3003 27027 v 1 99 v ; 27027 v ' 1287 J 

f 1312tt 4 . .... 119488tt 4 . 2 ., M f34887r 4 .. 2 224tt 4 . ... 82688tt 4 .A 

< uu u (uu) > u + < M H uu H u > u 

\ 2079 9009 v ; J \ 10395 495 135135 J 



d7? iol3 f 31047T 4 ., 256tt 4 . ...,o 9056tt 4 .. 4 8327r 4 .. 2 105152tt 4 .,. . 3485408tt 4 , 2 

-{uu) H tt u — — — — u {uu) — — — — u u + 



dr { 231 63 v ' 3465 945 3003 v ' 45045 

637568tt 4 . 4 ...... 6464tt 4 .. 2 ...... 86528tt 4 , . ...... ... . 35264tt 4 ,,2 30208tt 4 . 2 . ....... 512tt 4 . . . 

+ ~l5MT U {UU) - ^T" {UU) + -M6T {UU){UU) + ~M6T U " ~ ~M6T U {uu) -^^" )( ""» \ 1 

64tt 4 . ........ 2304tt 4 . ..... .... 29504tt 4 237376tt 4 . . ... 63040tt 4 . 2 ...... . 1088tt 4 .... . , 

(uu) (uu) (uu) u (uu) u (uu) H u (uu) H u (uu) > w M + 

33 y ' 143 v A ' 9009 v ; 9009 v ' 9009 v ' 117 

f359687r 4 60160tt 4 . .. n2 3968?r 4 ........ 58048tt 4 .,, . .... 121600tt 4 . 2 ,1 ..„ 

(uu ) ( uu ) H — — - — u (uu) — — — u 2 u 2 } u^- 



I 



3861 9009 v ' 2079 v ' 9009 v ' 9009 

f6407r 4 •• 134144tt 4 2 , 1 -fx f48647r 4 2 6592tt 4 • 4864tt 4 , 
— s uu — — — — — u (uu) > u + < u + uu H — u > u 



{ 693 9009 v ') \ 10395 10395 6435 



io (1 ° |0) = - ,f T M [82368k? + 80652«;f«;| + 12012/t 4 «;| + Q\4Qk\k%k\ + 143(k?k^ + k\k\k\k\ + 2k\k%k\ + k\i^k\)\ 

P o 10|1) = I!!! [51480k? + 97812k?k^ + 25311k 4 k 2 + 1387i4k2 k 3 + 572(k?k 2 ; + k\k\k\k\ + 1k\k%k\ + k\kIk\)] 
loo loo 

P (10|2) = ^ [-23634k? + 247104^^ + 99099k 4 k^ + 53911k 4 k^k^ + 2431(k?k 2 ; + + + k^k^)] 

p o 10l9) = tJtt^ [35100k? + 574860k?k^ + 210210k 4 k^ + 113971^^ + 5005(k?k^ + + 2k\k\k\ + )1 

135135 

p a m = ,or!"or [-302419k? + 108537k?k^ + 133419k 4 k^ + 72358^^ + 3718(k?k^ + k\k\k\k\ + 2k\k\k\ + k?k^k 4 )1 
135135 
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327T 4 

ho\o = ~ ,or,or (66524k?k 2 ! + 150140k?k? + 11336^^ - 2704k?k 2 k 2 - 22048k?kik^ + 18278kikik 2 - 8684k 2 kiKi+ 
135135 



+6903/«i + 5629^14 + 55872k? - 41904k?ki + 12207k 2 k 2 + 37440k?kik 2 k 2 + 2925k?k 2 + 48880k 2 k 2 k 2 

+ 143[(«i — 3ki«2 ~ 3«iK2 K 2) 2 + (3«iK2 + 3«iK 2 + Ki'K 2 — K1K2) 2 + (2k 1 K 2 K 3 + K1K2K3 + 3KiK 2 K3) 2 + 
+ 444(^4 + K 3 + K D — 2kiK 2 K 2 (3kiK 2 + K1K2 + 3kiK 2 )]) 



647T^ 

7io|i = _^ (82152k?K2 + 109712k?k 2 + 24076^^ ~ 5226k?k 2 k 2 - 41028k?kiK 2 + 15600kiKik 2 - 10816k 2 kiKi+ 

+5200/^ + 12701k^k^k^ + 36000k? - 32400k?ki + 11700k 2 k 2 + 64220k?kik 2 k 2 + 7475k^k 2 + 92248k 2 k 2 k 2 ,+ 

+572[(ki — 3kiK 2 — 3K1K2K2) 2 + (3kiK2 + 3kik 2 + kik 2 — K1K2) 2 + (2kik 2 K3 + kik 2 K3 + 3kik 2 K3) 2 + 
+k\k%k\{2k% + K 3 + K l) ~ 2kik 2 k 2 (3kik 2 + kik 2 + 3kik 2 )]) 



32tt 4 

J 10 , 2 = (187782k?*; 2 , + 8980k?k? + 93483k?/4 - 18460KfK 2 K 2 - 139828k?kik 2 ! + 13754kikik? - 31356k?kiki- 

1 135135 " 

+3250/^ + 48711k^k 2 .k 2 ; - 44606k? - 2478k?k x + 12675k 2 k 2 + 232544k?ki k 2 k 2 + 30251/tiK 2 , + 353522k 2 k 2 k 2 + 



+2431[(ki — 3/iiK2 — 3KiK2K2) 2 + (3kiK2 + 3kiK2 + K\K2 — KiKj) 2 + (2K1K2K3 + K1K2K3 + 3KiK2K 3 ) 2 + 
+ K 2 K 2 K 2 (2K2 + K 2 . + K 2 ) — 2kiK 2 K 2 (3kiK2 + KlK 2 + 3/ilK 2 )]) 



32tt 4 

J 10 |g = (221044k?k 2 + 84050k^k 2 + 99\h\n\n\ - 19812k*k 2 k 2 - 150852k?kik 2 + 22893kiKik 2 - 35698k 2 kiKi- 

_LoO loO 

13403. 4 103051 . o 2 o ... 37557 „„ „„, 3. . 63427 ,., „„„„„„ 2 . 2 2 
H — k^H — k^k 2 k 2 -16670k?-23430k?kiH — k 2 k 2 +251264k?kik 2 k 2 H — k^k 2 ! +377962k 2 k 2 k 2 .+ 

+5005[(ki — 3/iiK2 — 3kik 2 k 2 ) 2 + (3kik 2 + 3kik 2 + kik 2 — K1K2) 2 + (2kik 2 k 3 + kik 2 k 3 + 3kik 2 k 3 ) 2 + 

+ K 2 K2K 2 (2K2 + Kg + K 2 ) — 2kiK 2 K 2 (3kiK2 + KlK 2 + 3/tlK 2 )]) 
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327T 4 

(38659k?k^ - 677225k^ 2 + 124657/c^ - 21320k^k 2 k 2 - 151268kJki^ - 70616kiKik 2 - 20488k 2 kiKi- 



135135 



-28353/^ + 64298^^-263119^ + 182538^1-4438 

+3718[(ki — Skini, — 3K1K2K2) 2 + (3kiK2 + 3kiK2 + K1K2 — kiKj) 2 + (2K1K2K3 + K1K2K3 + 3/iiK2K3) 2 + 
+kIk2kI(2k2 + k| + «1) — 2k 1 k 2 k 2 (3k 1 K2 + KlK 2 + 3KiK 2 )]) 
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Appendix: MAPLE programs 

This Appendix contains a MAPLE program that allows one to obtain the quantities listed in the tables with 
the help of a computer. 

MAPLE program 

The MAPLE program we used in calculations looks as follows (note that uu[i,j] in the program is nothing but 
Uij, A\p, q] is A P ,U 2p ~ q+1 , is Ki(r) etc.; P[d, m]'s are the angular integrals, formula (16)): 

> #input I -> L then input P and apply L to P then input uu(k) and convert into k 
> 

> # u[i] = d~iu: uu[i,j] = (d~iu d~ju) : 
> 

> 

> M:=4: J:=3: 

> A:=array(0. .M,0. .M) : B:=array(0. .M,0. .M) : II : =array (0 . . J ,0 . . 2*M+2) : GG :=array(0. .M,0. . M) : 

> k:=array(0. . M) : K : =array (0 . . M) : N:=array(0. .M) : tu:=array(0. . M) : u:= array (0..M): 

> uu:=array(symmetric,0. .M,0. .M) : v:=array(0. .M) : VV:=array(0. .M,0. .M) : 

> for i from to M do: unassign( 'u [i] ' , 'uu[i,i] ' , 'tu[i] ') : 

> for j from to i-1 do: 

> unassign( 'uu [j , i] ' , 'VV[j,i]'): 

> od: od: 
> 

> # Parameters A[p,q] 

> W[0,0] :=1/V(y): 

> for p from to M-l do: 

> W[p+1,0] :=diff (W[p,0] ,y)/V(y) : 

> for q from to p do: 

> W[p+l,q+l]: = (diff (W[p,q+l],y) + W[p,q] )/V(y) : 

> od: od: 

> V(y) : =add(v [n] *y~n/n ! ,n=0 . .M) : 

> for p from to M do: for q from to M do: 

> A[p,q] :=simplify(v[0] ~(2*p+l-q)*eval(VV[p,q] ,y=0)) : 

> od:od: 
> 

> # list entries of TABLE II 

> for p to M do 

> for j from to p do 

> print (A [p, j] ) : 

> od: od: 
> 

> 

> # entries of TABLE III: Polynomials P[d,m] 

> 

> x:=t*t + a[0] *u+add(a[ll] *u[ll] ,11=1 . .M) : xu:= t*tu[0]+add(a[12] *uu[0,12] ,12=0. .M) : 

> x2:=t~2*t2 + a[0]~2 + add(a[nl] *a[nl] *uu[nl ,nl] ,nl=l . .M)+add(2*t*a[13] *tu [13] , 13=0 . .M)+ 

> 2*add(add(a[n]*a[14]*uu[n, 14] ,n=0. .14-1) ,14=1. .M) : xL:= t*Lt+a[0] *Lu+add(a[15] *Lu[15] ,15=1 
> 

> x; xu; x2; xL; # print for illustration purposes 

> 

> 

> for p to M do 

> for 1 to p do 

> P[2*p+2,2*1] :=2-(p+l)*Pi~p*2-(p-l+l)*(p-l+l) !/(2*p-l+2*l) ! *add( (-1) ~n*(2*l) !* 

> (xu)~ (2*l-2*n)*(x2)~n/2~n/n!/(2*l-2*n) !*2"(p-l+2*l-n)*(p-l+2*l-n) !/2~(p-l+l)/(p-l+l) ! ,n=0. 

> od: od: 
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> 

> for p to M do 

> for 1 from to p do 

> P[2*p+2, 2*1+1] :=2~(p+l)*Pi~p*2~(p-l+l)*(p-l+l) !/(2*p-l+2*l) ! *add( (-1) ~n* (2*1+1) !* 

> (xu)-(2*l+l-2*n)*(x2)~n/2~n/n!/(2*l+l-2*n) ! *2~ (1-n) * (p+2*l-n) !/(p+l) ! ,n=0. .1) : 

> od: od: 
> 

> # G[d,q,r] 

> 

> # expect long calculation 

> for p to M do 

> for q from to p do 

> for r from to p do 

> J J : =A [p , q] *A [p , r] *P [2*p+2 , 2*p+l-q-r] : 

> for j from to M do 

> LLL:=add(ii!*coeff (coeff (JJ,v[j] ,ii) ,a[j] ,ii) ,ii=0. .2*M) : 

> JJ:=LLL: 

> od: 

> G[2*p+2,q,r] :=diff (LLL,t)/(2*p+l-q-r) ! : 

> od: od: 

> od: 



> # expect long calculation 

> 

> # relations between uu 

> 

> uu[0,0]:=l: 

> for j to M do 

> uuu:=add(j !*uu[n, j-n] /n! / (j-n) ! ,n=l. . j-1) :uu[0, j] := -uuu/2: 

> uu[0,j] := -uuu/2: 

> od: 
> 

> # entries of TABLE V: intensities of radiation of spin-J particles 

> 



> for p to M do 

> for q from to p do 

> for r from to p do 

> GG[q,r] :=simplify(G[2*p+2,q,r] ) : #print(simplify(G[2*p+2,q,r] ) ,2*p+2,q,r,GG[q,r] ) : 

> od: od: 

> II[0,2*p+2] :=GG[0,0] : 



> for j to J do 

> for q from to p do 

> for r from to p do 

> GG[q,r]:= add(add( uu [ql ,rl] *GG [q+ql ,r+rl] ,ql=0 . .p-q) ,rl=0 . .p-r) : 

> unassign( 'ql ' , 'rl ' ) : 

> od: od: 

> II [j ,2*p+2] :=simplify(GG[0,0] ) : 

> od: 



> od: 
> 

> for p to M do 

> for j from to J do 

> print (2*p+2, j ,11 [j ,2*p+2] ) : 

> od: od: 
> 

> # FRENET CALCULATION (TABLE IV) 

> # expect long calculation 

> uu[0,0] :=l:uu[l,0] :=0:uu[l,l] :=-k[l] (z)~2: N [0] :=u[0] :N[1] :=l/k[l] (z)*u[l] : 

> k[0](z):=l: K[0](z):=l:K[l](z):=k[l](z): 
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> for p from 1 to M-l do: 

> uu[p,p+l]: = diff (uu[p,p] ,z)/2: 



> K[p+1] (z) :=-K[p] (z)*k[p+l] (z) : 

> for i from to p-1 do: 

> uu[i,p+l] :=-mi[i+l,p]+diff (uu[i,p] ,z) 

> od: 

> B[p+1,0] :=uu[0,p+l] : 

> for i from 1 to p do 

> B[p+l,i] :=-add(uu[jj ,p+l] *dif f (N [i] ,u[jj] ) , jj=0. .p) : 

> od: 

> N[p+1] : = (u[p+l]-add(B[p+l,ii]*N[ii] ,ii=0. .p))/K[p+l] (z) : 

> uu[p+l,p+l] := uu[0,p+l] ~2 - K[p+l](z)~2 -add(B [p+1 , ii] ~2, ii=l . .p) : od: 

> # expect long calculation 

> 

> # list entries of TABLE IV 



> for p from to M do 

> for j from to p do imu[p, j] :=simplify(uu[p, j] ) : 

> print (mm [p, j] ) : 

> od: NN[p] :=simplify(N[p]) : print (p,NN[p]) : od: 

> # end of FRENET 

> 

> # intensities through FRENET (TABLE V) 

> for p to M do 

> for j from to J do 

> III [j ,2*p+2] :=simplify(II [j ,2*p+2]) : 

> print (2*p+2,j ,111 [j,2*p+2]) : 

> od: od: 

> # P_0 — projection to u (TABLE V) 

> for p to M do 

> for j from to J do 

> Hp [ j , 2*p+2] : =subs (tu [1] =0 , tu [2] =0 , tu [3] =0 , tu [4] =0 , III [j , 2*p+2] ) : 

> print (2*p+2, j , Hp [j ,2*p+2] ) : 

> od: od: 

> # contraction with u: t=u, tu=uu (TABLE V) 

> 

> for j from to M do 

> tu[j] :=uu[0, j] : 

> od: 

> for p to M do 

> for j from to J do 

> IIIILj ,2*p+2] :=simplify(II[j ,2*p+2]) : 

> print (2*p+2, j ,1111 [j ,2*p+2] ) : 

> od: od: 

> for j from to M do 

> unassign( 'tu[j] ' ) : 

> od: 

> # end 



Prenet basis 

A few comments are needed to explain how the Frenet curvatures are handled with the program. 

We deal with the Frenet basis recursively. Our main objects are the scalars Uij, the coefficients Dij of expand- 
ing Frenet basis, iVW = J^_ Dikd^u and the coefficients Bij of the inverse expansion d° T u = X)i=o A?',*-^"^- 
All these quantities are dealt with as functions of the Frenet curvatures and their derivatives. 

Now, the knowledge of these quantities for alH, j < p allows one to calculate them for alH, j < p+ 1 in the 
following way: 
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Calculate B p+ i p+ i = —K p+ iB pp . To prove this formula, one suffices to note that d T N^ = -^-dP +1 u + 

. . . , which has to be expanded into KpN^" 1 ^ — n p +iN^ p+1 ^ . However, d p+1 u contributes to only. 
Therefore, d T N^ = B "+ 1 '"+ 1 AKP+ 1 ) + . . . and k v+1 = - Bp + 1 - P+1 . 

PiP P>P 

Calculate B iyP+ i for all < i < p — 1 using 

Ui, p +i — d T Ui iP — Ui + i. p 

and U p . p +i = \d T U pp . 

Calculate Bi tP+ i for all 1 < i < p using Bi tP+ i = —N^dv +1 u = — ^2 k DikUk, P +i, the latter expression 
can be simbolically written as j— J2 k q^I^ t^fc,p+i|- Note that -B p +i,o = Uo tP +i- 

Calculate = y=° P+ " . 

Calculate U p+lyP+1 = (B p+1 . ) 2 - X)f=i (#p+m) 2 - 
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